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We study power corrections to the event shapes in e + e~ —annihilation in the two jets kinemat- 
ical region, e<l. We argue that for e ~ Aqcd/Q all power corrections of the form l/(Qe)" 
have to be taken into account for an arbitrary n. This is achieved by introducing a new 
universal distribution, the shape function, which describes the energy flow in the final state. 
The event shape differential distributions are given by a convolution of the shape function 
with perturbatively resummed cross-sections. Choosing a physically motivated ansatz for the 
shape function, we observe a good agreement of our predictions with available data on the 
differential thrust, C-parameter and heavy-jet mass distributions and their first few moments. 



Introduction 



The experimental data for the event shape distributions in e + e~— annihilation deviate from 
perturbative QCD predictions by corrections suppressed by a power of the hard scale l/Q p 
ON ■ with the exponent p depending on the observable. These power corrections are associated 



with hadronization effects in the final states. Their contribution to the differential event shape 
distributions and the mean values has a different form. The leading hadronization corrections 
to the mean values can be parameterized by a single nonperturbative scale, whereas for the 
differential distribution in the end-point region, e ~ A-qcd/Q, one has to take into account 
an infinite set of the power corrections of the form \/{eQ) n for arbitrary n. Denoting by e a 
generic event shape variable (e = 1 — T, p, C), one can write the expression for the differential 
distribution with perturbative and hadronization corrections included in the following general 



_l_da = da^ + l Pwr(fte) . 
a tot de de Q p 
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Here, Ff mi j -r (Q,e) receives both perturbative and nonperturbative contributions, which need to 
be disentangled in order to understand the physical origin of power corrections. We shall study 
the event shape distributions in the two-jet region < e < e max , with the upper limit e max 
separating the regions with two and three jets in the final states. For the event shapes under 
consideration pQCD predictions are known to two-loop accuracy: 

^f = ^A e {e)0(e max -e) + (^)^e(e) + O(a»(0)) (1) 

with A e (e) and B e (e) being coefficient functions. Due to enhancement of soft and collinear gluon 
emissions, these functions become divergent in the end point region, e — ► 0, 

with to = po = 1 and Co = 6. Moreover, large Sudakov logs 0^ In 2 ~ n (e)/e, n > are present 
to all orders of perturbation theory. They spoil the convergence ofjperturbative expressions 
and need to be resummed. This has been done to the NLL accuracy □. It turned out that the 
resummed pQCD predictions do not fit the data in the end-point region even at the Z-mass scale 
(Fig. |l]a), thus calling for a better understanding of the underlying nonperturbative physics. 
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2 Event shape distributions 

To understand the structure of the power corrections to the event shape distributions, we 

first have to identify the relevant scales. In the end-point region, the final states in the 

e + e~— annihilation consist of two narrow quark jets surrounding by soft gluon radiation. Going 

through QCD analysis of such states, one finds the dynamics is driven by two different scales, 

the so-called soft scale, Qe, and collinear one, Q\/e. The soft scale Qe sets up a typical momenta 

of the soft gluons, while the collinear scale Q\fe determines the transverse size of the outgoing 

fast jets. Nonperturbative corrections to the distributions appear suppressed by powers of both 

scales. We notice, however, that in the end point region, e ~ Q q D , the following hierarchy of 

the scales holds, Qe <C Q\/e <C Q. This suggests to neglect power corrections on the larger scale 

Qy/e and resum all corrections on the smallest scale Qe. The resulting differential distribution 

can be written as 

1 da , ,„ x , , „ „ ,„ „ ,„„ , 2 



a (a s (Q),ln(e),l/Qe)+O(l/Q z e). (2) 

&tot de 

This approximation amounts to neglecting the internal structure of two jets and treating them as 
two fast classical color charges moving in the direction of jet momenta and emitting soft gluons in 
the final states. The resulting expression for the differential distribution, Co(a s (Q), ln(e), 1/Qe), 
can be described using Wilson loop approach UU. By the construction, it should resum both 
large perturbative Sudakov logarithms and power corrections on the smallest scale Qe 

°o = ^ + tj^Z k (a s (Q)M(e)). (3) 

Here, A& are some nonperturbative scales depending on the choice of the event shape variable, 
and Efc(a s (Q),ln(e)) are perturbatively calculable coefficient functions. The equation (||) can 
be thought of as a generalization of the OPE expansion to the event shape distributions. For 
e ^> Xqcd/Q the leading nonperturbative corrections come only from the first term in the sum 
(IH). Their net effect on the distribution is the shift of the pQCD spectrum towards larger values 
of the event shape trl 



v 



o"o= j^-(e-Ai/Q), as e> Xqcd/Q- 

a to t de 

The IR renormalons models give the same prediction with Ai being universal nonperturbative 
parameter related to the integral of the coupling constant over small momentum region. 

For e ~ Q q D all terms in (|3|) become equally important and need to be resummed. This 
is achieved by introducing a new nonperturbative distribution, the so-called shape function. 
Taking the thrust as an example we have |J 

1 drr rQt drr pert 

~ f def t (e)^-(t-e/Q). (4) 



crtot dt Jo dt 

Here, /t(e) = / deaden /(e^, en)5(e — el — en) with the shape function /(e^, e^) being a universal 
nonperturbative distribution describing the energy flow into the left and right hemispheres in 
e + e~— annihilation final states. The universality of /(ez,, er) implies that it does not depend on 
the hard scale Q. 

The shape function receives contributions of two different types 

/(ei, e R ) = /incl( e ^)/incl( e R) + ^non-incl^' e «) ( 5 ) 



1 Similar formulae hold for the other event shape variables. 



The first one is the inclusive contribution from gluons produced and decaying into the same 
hemisphere. This kind of event is partially taken into account by IR renormalons models. The 
second type of contribution corresponds to non-inclusive events, that is when gluon is produced 
in one hemisphere but decays into particles flowing into another hemisphere. The non-inclusive 
contribution describes the cross talk between two hemispheres and we expect that this effect is 
important for non-inclusive variables like the heavy-jet mass. 

The shape function is a well defined nonperturbative QCD distribution Era which in present 
can not be calculated from the first principles. Therefore, confronting our predictions with the 
data we shall rely on particular ansatz for this function. Choosing the ansatz, we require that 
the shape function should be positively definite, vanish at high energies and have a power like 
behaviour at the origin due to the phase space suppression. This leads to 
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M[a,b) (e L e R y 
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exp 
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(6) 



Here, M(a, b) is the normalization constant, a, b and A are nonperturbative parameters that 
have to be fixed by comparing our predictions for the event shape distributions with the data. 
The parameter a determines how fast the shape function vanishes at the origin, A sets up a 
typical scale of the soft radiation and b controls the size of non-inclusive corrections, so that 
^/non-incl( eL ' £r ) = at & = 0. We fix the parameters by confronting the predictions for the 
heavy-jet-mass and C— parameter distributions with the data at Q = Mz (see Fig. [l], where we 
present the predictions for heavy-jet mass distribution.) 
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A = 0.55 GeV. 



To test universality of the shape function, one can use these valuesjjf the parameters to compare 
the same distributions with the data over a wide energy interval EJ, 35 GeV < Q < 189 GeV. 
As was shown intJ'El, a good agreement is observed over the whole range of the event shapes 
including the end-point region e ~ Aqcd/Q- 
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Figure 1: (a) Heavy-jet mass distributions at Q — Mz with and without power corrections included, 
(b) Comparaison of the QCD predictions for heavy jet mass at different energies from top to bottom Q(GeV) = 

35, 44, 91, 133, 161, 172, 183. 
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Figure 2: Comparison of the QCD predictions for the second moments (p ) and (C ) with the data. Dotted 
lines denote 0(a 2 s )— perturbative contribution, solid lines take into account power corrections given by Eq. bj 



3 Moments 



Using the obtained expressions for the differential distribution for the heavy-jet mass and the C- 
parameter, we can calculate their lowest moments. This calculation can be performed provided 
that the moments are dominated by the contribution of two-jet final states. One can verify that 
this is the case for the first and the second moments El. The first two moments of various event 
shape distributions have been measured by the LEP experimentsu. The obtained values deviate 
from the predictions of IR renormalon models, especially for the second moments. 

Going through the calculation of the second moment of the C-parameter and heavy-jet mass 
distribution we find 
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where A 2 is equal to the second moment of the shape function, 
measures non-inclusive contribution to the heavy-jet mass, 



(e 2 ) = 1.7 GeV 2 and <5A 2 



5X 2 



<(cl - crY 



■pmax f dapr\ 2 

1+4 / dpp 1 

o 



\ dp ) 



Here, ((e L - e R ) 2 ) = 0.14 GeV 2 and ^ max dpp (j^-) varies from 2.19 for Q = 10 GeV to 1.85 
for Q = 100 GeV. We notice that non-inclusive correction <5A 2 does not affect the second moment 
of the C— parameter. The reason for this is that the heavy-jet mass is less inclusive observable 
with respect to a single jet than the C— parameter and, therefore, it is more sensitive to the 
cross-talk effects between two hemispheres. Non-inclusive effects provide important contribution 
to the second moment of the heavy-jet mass ~ 14%. These corrections have not been taken into 
account in the IR-renormalons models, since there the two hemispheres were supposed to be 
uncorrelated. Finally, using (]?]) we observe a good agreement of the QCD predictions with the 
available data for the second moments of the heavy-jet and C— parameter distributions (see 
Fig- (!))• 



4 Conclusion 

The shape function approach allows to reveal the physical meaning of the nonperturbative power 
corrections to the event shape distributions in the e + e~— final states. The emerging structure 
of power corrections to the event shapes looks as follows. For e 3> Aqcd/Q the main nonper- 
turbative effect is the shift of the pQCD spectrum towards larger values of the shape variable, 
e — ► e — A e /Q, with the scale A e given by the first moment of the shape function. In the end 
point region, e ~ Aqcd/Q, all power corrections ~ l/(Qe) n become equally important. They 
can be resummed in a new distribution - the shape function which is a universal nonperturba- 
tive distribution describing the energy flow into two hemispheres in the final state. It takes into 
account both inclusive and non-inclusive corrections with the latter being especially important 
for the observable like heavy-jet mass, which are not completely inclusive with respect to a single 
jet. Choosing ansatz for the shape function, we observed a good agreement with the data for the 
differential thrust, C— parameter and heavy-jet mass distributions as well as their mean values 
and the second moments. 
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